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We perform ab initio calculations for nuclei around 16 0, 40 Ca, and 56 Ni using realistic nucleon- 
nucleon forces. In particular, Ni is computed as the heaviest nucleus in this kind of ab initio 
calculation. Ground-state and single-particle energies including three-body-cluster effects are ob- 
tained within the framework of the unitary-model-operator approach. It is shown that the CD-Bonn 
nucleon-nucleon potential gives quite good results close to the experimental values for all nuclei in 
the present work. 
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One of the most fundamental problems in nuclear 
physics is to describe and understand nuclear proper- 
ties from the underlying nuclear forces. To solve this 
problem, a realistic nucleon-nucleon (AW) interaction [l[, 
which has a strong repulsive core and a complicated spin- 
isospin structure, has been employed Q. In addition, to 
obtain morequantitative results, a three-nucleon (NNN) 
interaction [J [3] has been used in some cases. Light nu- 
clei with the mass number up to A ~ 12 have been well 
understood from the ab initio calculations employing the 
realistic NN and NNN interactions with the Green's func- 
tion Monte Carlo (GFMC) method [111] and tne no-core 
shell model (NCSM) d, Q. While these methods have 
been successful in the light nuclei, their applications to 
heavier systems become rather difficult due to exponen- 
tial increase of the computer performance to be needed. 

Coupled-cluster (CC) theory or, in other words, e s (or 
e T ) method HQ is a promising one for the microscopic 
calculation for the heavier nuclei. Recently, the first cal- 
culations for nuclei up to the p/-shell region, 48 Ca and 
48 Ni, have been reported with a CC method including 
the excitations of the singles and doubles which is re- 
ferred to as CCSD [lfjj|. Well-converged results in a suffi- 
ciently large model space have been obtained using a chi- 
ral N 3 LO NN interaction [ll[ as one of the realistic NN 
forces. While the CCSD calculations give results fairly 
close to the experimental values, there still remain some 
discrepancies between the results and experiments. One 
of the reasons of the discrepancies may be attributed to 
the missing AW/Vinteraction in the calculation. Although 
the NNN force has been considered to be an indispensable 
ingredient for a more quantitative description of the nu- 
clear properties, there has been no definite way of using 
the NNN force directly in the calculation for the heavier 
nuclei. 

Given this situation, it is still worthy to compute nu- 
clear properties using only the realistic NN interaction 
in a rigorous way and to investigate to what extent nu- 
clei can be described with only the AW force. Such a 
study could be helpful to evaluate the magnitude of the 



NNN-ioice effect in heavier nuclei in future works. 

In this Letter, we report the results of calculated 
ground-state energies and single-particle ones for hole 
states in nuclei around 16 0, 40 Ca, and 56 Ni with the 
unitary-model-operator approach (UMOA) [l~2, HI, 14 1. 
The calculation for 56 Ni, which is a typical p/-shell nu- 
cleus, is performed for the first time within the UMOA 
framework. Furthermore, 56 Ni is the heaviest nucleus 
for which this kind of ab initio calculation has been per- 
formed. In the UMOA, a Hermitian effective interaction 
is derived through a unitary transformation [15l . Ilq ] . The 
unitary transformation method has been widely used in 
other microscopic methods in nuclear physics, such as the 
NCSM 0, 0] and the hyperspherical harmonics effective 
interaction method (EIHH) [17| . The unitary transfor- 
mation treats successfully short-range correlations due to 
the strong repulsive core of the AW force in a truncated 
model space, but the model space should be sufficiently 
large in the sense of the ab initio calculation. 

In the UMOA, a unitarily transformed Hamiltonian H 
of the original many-body Hamiltonian H is given in a 
cluster-expansion form as H = e~ s He s = + H^ 1 + 
H(3) + where S is a two-body anti-Hermitian oper- 
ator and is determined by solving a decoupling equation 
between the model space and its compliment [l8|. The 
Jji 1 ) 7 H^ 2 \ and H^ 3 ' are the one-, two-, and three-body 
cluster (3BC) terms, respectively. The method of the 
actual calculation and the results of nuclei around 16 O 
including up to the two-body cluster terms using mod- 
ern AW forces have been given in detail in our previous 
study [14]. In the present work, we apply this method 
to the heavier nuclei up to 56 Ni and evaluate effects of 
the 3BC terms systematically. As for the evaluation 
of the 3BC terms, we follow the prescription given in 
Refs. dlGl. 

In Fig.[TJ we first demonstrate HQ and p\ dependences 
of the calculated ground-state energies of 16 O includ- 
ing the 3BC effects. Here, h£l is the harmonic-oscillator 
(h.o.) energy of the single-particle basis states. The pi 
stands for a boundary number defined with a set of the 
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FIG. 1: (Color online) The hil and p\ dependences of the 
calculated ground-state energies Sg.s. of 16 0. 



h.o. quantum numbers \n\,l\\ and {n2,h} of the two- 
body states as p\ — 2ni + l\ + 2n 2 + h, and specifies the 
size of the model space of the two-body states. The Nijm- 

and the chiral N 3 LO 



I [191, the CD-Bonn [20j, and the chiral N 3 LO [ll| NN 
forces are employed as the realistic AW interactions. The 
Coulomb force is added to the proton-proton channel. It 
is seen that well-converged results with respect to the 
boundary number p\ are obtained at p\ — 18 and 14 for 
the Nijm-I and the CD-Bonn interactions, respectively. 
For the N 3 LO interaction, the convergence property is 
rather different from the other two forces. For example, 
at Ml = 15 MeV, it is difficult to search for convergence 
up to p\ = 14. However, for the larger values of pi, the 
results rapidly converge toward the point at p x =20. 

In TablelH we tabulate the energies of the one- and two- 
body cluster terms B( 1 + 2BC ), the 3BC terms E^ 3BC \ and 
the total ground-state energy E g _ s . of 16 0. The binding- 
energy per nucleon BE/A = —E g . a JA is also given. The 
values of Ml = 14 MeV and p x = 18 for Nijm I, Ml = 15 
MeV and p x = 14 for CD Bonn, and Ml = 15 MeV 



TABLE I: The calculated energies of the one- and two-body 
cluster terms £ (1+2BC) , the 3BC terms £( 3BC > and the total 



ground-state energ; 
taken from Ref. [21 



-Eg.s. of O. The experimental value is 
. All energies are in MeV. 
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FIG. 2: (Color online) The TiQ and pi dependences of the 
calculated energies E(= —BE) of the lowest 1/2" and 3/2~ 
states of the spin-orbit doublet in 15 O. 



and pi — 20 for N 3 LO are shown as the optimal ones. 
It is seen that, although the 3BC terms have attractive 
and sizable contributions to the ground-state energy, the 
calculated ground-state energies are still less bound than 
the experimental value. In the present calculation, a gen- 
uine NNN force is not taken into account. The inclusion 
of the NNN force could compensate for the discrepancies 
between the theoretical and experimental values, which 
has been shown in the microscopic studies of light nu- 
clei 0, [22|]. Note, however, that the energies of 93.6 % 
to the experimental value are attained from only the NN 
force for the CD-Bonn and the N 3 LO potentials, though 
84.4 % for the Nijm-I interaction. 

The 3BC effect for N 3 LO is significantly larger than 
the ones for Nijm I and CD Bonn. A similar tendency 
is seen in the recent ACCSD(T) computation including 



TABLE II: The calculated energies of the one- and two-body 
cluster terms _B( 1+2BC ), the 3BC terms £ (3BC) , and the total 
energy E(= —BE) of the lowest l/2~ and 3/2~ states of 
the spin-orbit doublets in ls O and 15 N. The quantity -E s . . is 
the spin-orbit splitting energy including the 3BC effects. The 
energy difference -Bdiff of the ground-state energies between 
15 O and 15 N is also given. All energies are in MeV. 
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FIG. 3: (Color online) Same as Fig. Q] but for 40 Ca. 



triples corrections [23| . We have found that the large 
3BC contribution also applies to the other nuclei in the 
present study. Owing to this property, we have not yet 
obtained the converged results for the other nuclei. For 
this reason, we do not show the other results for N 3 LO. 
One may also notice that the result of E^ 1+2BC ^ for N 3 LO 
shows a large difference of about 4 MeV from that in our 
previous study [l4[. This is due to strong dependences 
of J B( 1 + 2BC ) and £( 3BC ) on Ml and Pl for N 3 LO. 

In Fig. [2 we illustrate the Ml and p\ dependences of 
the total energy E(= —BE) including the 3BC effects of 
the lowest l/2~ and 3/2~ states of the spin-orbit doublet 
in 15 O. These states are representative single- hole states 
of neutron in 15 0. For Nijm I, we take the values of hfl = 

14 and 13 MeV at p x = 18 for the l/2~ and 3/2~ states, 
respectively, as the optimal ones, and for CD Bonn, Ml = 

15 and 14 MeV at p\ = 14. These optimal values are 
tabulated in Table [ill The results for the proton-hole 
states in 15 N are also given. The convergence properties 
for 15 N are similar to the case of 15 0. 

The microscopic description of the spin-orbit splitting 
in nuclei is a long-standing problem. In Table [TTJ the 
spin-orbit splitting is denoted by £7 S . . which is the dif- 
ference of the binding energies between the 1/2 - and 
3/2~ states. Our results show smaller splitting energies 
than the experimental values for these hole states, which 
does not contradict previous studies @, HI . We should 
note, however, that the magnitude of the lack of the split- 
ting energy depends considerably on the interactions em- 
ployed. For CD Bonn, the differences are only 0.34 and 
0.48 MeV for 15 and 15 N, respectively. The calcula- 
tion including the NNN force could give a better result 
as shown in Refs. 24, 25|. 
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FIG. 4: (a) (Color online) Same as Fig. [TJ but for 56 Ni. (b) 
The extrapolation curve of E g . B . of 56 Ni. See text for details. 



The energy difference of the ground states between 15 O 
and 15 N is denoted by -Ediff in Table HT1 Our results are in 
good agreement with the experiment. Similar tendency 
has been found in the case of 3 He and 3 H in our previ- 
ous work 



14] 



Since we include the Coulomb force, the 
small differences between the results and experiment may 
be attributed to the effects of the charge-independence 
breaking of the original NN forces. 

In Fig. [3l we show the Ml and p\ dependences of the 
calculated ground-state energies including the 3BC ef- 
fects of 40 Ca. We take the values of Ml = 13 MeV and 
pi = 20 for Nijm I, and Ml = 14 MeV and p x = 18 for 
CD Bonn as the optimal values. Since we handle a heav- 
ier system 40 Ca than 16 0, we need a larger model space 
to obtain the converged results. The optimal values are 



TABLE III: Same as Table U but for 40 Ca. 



given in Table IIIII The results of E g . s . are less attractive 
than the experimental value similarly to the case of 16 O. 
However, for CD Bonn, the calculation attains 99.5 % of 
the experimental energy, and the difference between the 
result and experiment is only 1.78 MeV. This difference 
is much smaller than that for 16 O despite the fact that 
the absolute value of the ground-state energy of 40 Ca is 
much larger than that of 16 O. This fact suggests that the 
NNN force plays a complicated role in the inner (dense) 
and outer (thin) regions of the nuclei. 

In Fig. SJa) , the Ml and p\ dependences of the calcu- 
lated ground-state energies including the 3BC effects of 
56 Ni are illustrated for CD Bonn. A converging result is 
seen at Ml = 14 MeV and pi = 18. However, the en- 
ergy difference between the two values for pi = 16 and 



TABLE IV: Same as TableQ] but for 56 Ni. The results only for 
the CD-Bonn potential are shown. The extrapolated values 
for pi — » oo are also given. 
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18 at the energy minima amounts to 7.38 MeV which is 
somewhat large compared to the case of 40 Ca where the 
difference is 1.97 MeV. The results for the CD-Bonn po- 
tential show a regular pattern of convergence (in contrast 
to the results for N 3 LO in Fig. |T|) . In order to estimate 
the remaining effect of the larger model space, we per- 
form an extrapolation, as given in Refs. [jj, using 
the following formula: E g _ s _(pi) = -Boo + Ce~ 7p i, where 
.Boo j Cj an d 7 are the coefficients determined in the least- 
squares fitting procedure. We have found that the data 
points for hfl — 14 MeV from p\ = 8 to 18 are well fit- 
ted with this formula. The curve given by the formula 
is shown in Fig. |4jb) . The optimal values of the coeffi- 
cients are £«, = -473.17 MeV, C = -316.06 MeV, and 
7 = 1.3148 x 10~ 2 . Therefore, the extrapolated ground- 
state energy for p\ — ► oo becomes E g . a ,(pi —> oo) = 
Eoo = —473.17 MeV. The difference between the extrap- 
olated value and the result for p\ — 18 is 3.82 MeV, and 
thus the result for p\ = 18 is considered to be an almost 
converged value. The extrapolated value reproduces 97.8 
% of the experimental ground-state energy. 

In Table 033 the optimal value of TiQ, = 14 MeV and 
pi = 18 and the extrapolated one for 56 Ni using the 
CD-Bonn potential are listed. The extrapolation for 

£(1+2BC) m the 

same manner has been performed, and 
its result is also given. It is seen that the 3BC effect of 
56 Ni is considerably larger than that of 40 Ca. This may 
reflect the difference of the shell closure, namely, O/7/2 
sub-shell closed for 56 Ni and IsOd major-shell closed for 
40 Ca. 

In summary, we have applied the UMOA to the ground 
states of the closed-shell nuclei 16 0, 40 Ca, and 56 Ni, and 
the single-hole states in 15 O and 15 N. The p/-shell nu- 
cleus 56 Ni is the heaviest one for which this kind of ab 
initio calculation has been performed. The binding ener- 
gies including the 3BC effects have been obtained using 
the Nijm-I, the CD-Bonn, and the chiral N 3 LO AW inter- 
actions. We have found that the chiral N 3 LO interaction 
gives rather large 3BC contribution to the ground-state 
energy of 16 O compared to the other two forces. All re- 
sults lack the binding energies in reproducing the experi- 
mental data. However, the magnitude of the missing en- 
ergy depends considerably on the interactions employed. 
The CD-Bonn potential gives quite good results close to 
the experimental values for all nuclei in the present work. 
The inclusion of the NNN force is expected to make at- 
tractive contributions and compensate for the remaining 
discrepancies between the results and experiments. 

One may compare t he p resent results with recent mi- 
croscopic calculations [lfl IS IS 



fundamental forces is of great interest. Recently, a novel 
NN force from a lattice QCD calculation has been re- 
ported, and a more elaborate work of nuclear force is in 
progress [2^|. We will pursue studies using forces based 
upon the lattice QCD. 
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using realistic 

NN forces including the ones used here. Although the 
present and recent methods give similar results, there 
still remain some discrepancies. It is an important prob- 
lem to clarify the origins of the discrepancies in order to 
develop the microscopic many-body methods. 

For a deeper understanding of nuclei, the use of more 
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